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Abstract: Moments of secular and inverse secular coefficients, averaged over random 
matrices from classical groups, are related to the enumeration of non-negative ma¬ 
trices with prescribed row and column sums. Similar random matrix averages are 
related to certain configurations of vicious random walkers and to the enumeration 
of plane partitions. The combinatorial meaning of the average of the characteristic 
polynomial of random Hermitian and Wishart matrices is also investigated, and 
consequently several simple universality results are derived. 


1 Introduction 

The richness of random matrix theory was greatly enhanced during the last couple of years 
of the millennium by the discovery of its intimate connections with increasing subsequences 
and non-intersecting lattice paths in enumerative and asymptotic combinatorics (for reviews 
see QEC2) These topics relate to non-negative integer matrices via the celebrated Robinson- 
Schensted-Knuth (RSK) correspondence mum. A special class of non-negative matrices are 
so called magic squares — they have the property that the sum of the elements in any row or 
column is equal to a prescribed positive integer j say. Such matrices are natural objects in RSK 
theory, so one might expect a relationship between magic squares and random matrices. This 
is indeed the case, and its development forms the main theme of the recent work [3] (see JJJ for 
related results). One of the goals of the present paper is to continue the development of this 
theme. 

Magic squares are special cases of non-negative integer m x n rectangular matrices in which 
the sum of the elements in each row j is equal to fiji while the sum of the elements in each 
column k is equal to pk • Without loss of generality, for the purpose of enumeration we can 
insist that Hi > ■ ■ ■ > Hm and /A > • • • > p n so that /j, := (jii, • • •, Hm) and p := (pi, ..., p n ) 
form partitions. We denote the total number of such matrices by AV. Writing the partitions 
in terms of the frequencies of their parts by 

H = {l ai ---l ai ), p = (l hl ■■■l bl ) (1.1) 


1 


where l = max(/ti, n\), it was proved in 0, Thm. 2] that for N > max(^ j=] ja-j. Y2j=\j^j) 

i 

Nrf = E MeC 7(iV) n (S Cj (M)P(S Cj (M)) bj . (1.2) 

3 = 1 

In (II .21) the average is over matrices M chosen from the group U(N ) at random with respect 
to the Haar measure (uniform distribution), while Sc j(M) is the jth secular coefficient of the 
characteristic polynomial of M, 

N 

P M (z) = det(M - zI N ) = ^2 Sc j{M)(-z) N ~ j . (1.3) 

3=0 

From it follows that the number Flfc(j) of k x k magic squares, specified as k x k non¬ 
negative integer matrices with the sum of elements in each row and column equal to j, is for 
N > kj, given by the formula 

H k (j) =K M£U{ N)\Sc j (M)\ 2k . (1.4) 

We will extend the formula on in a number of directions. The first relates to counting 
formulas analogous to o in the case that the matrices are constrained by a symmetry property. 
Two such formulas, both relating to n X n symmetric matrices, are known from 0. Thus with 
the row sums (which must be equal to the column sums as the matrices are symmetric) labelled 
by the partition /r = (/j, i, ..., /j , n ) = (l ai ■ ■ - l ai ) we have that for \fj,\ := Y^j=i J a j even, N > |/r|, 
and all elements on the diagonal zero, the total number N® of such non-negative integer matrices 
is given by 

l 

< = ^MeO(N) II(Sc,(M)r. (1.5) 

3= 1 

If instead all elements on the diagonal are permitted to be even, the total number A^ p of such 
matrices is given by 

i 

= E MeC/Sp (27V) n(Sffi(M)r, (1.6) 

3=1 

In (11.51) the average is over the group of orthogonal matrices O(N), while in (11.61) it is over the 
group of 2N x 2N unitary symplectic matrices USp(2N) (both with the Haar measure). In 
Section 2 we will extend CC3 and CD to the case that the diagonal elements have no special 
restriction. We also give analogous formulas in the case of square matrices symmetric about both 
the diagonal and anti-diagonal, and for matrices symmetric about the centre point of the matrix. 
That such symmetrizations relate to averages over classical groups is a consequence of results of 
Rains m, and Baik and Rains |21 • Another tractable case in terms of a random matrix average 
to be considered is the setting relating to a but with the entries of the matrix restricted 
to 0’s and l’s. If the matrix has a 2 x 2 block structure, with the two diagonal blocks having 
non-negative integer entries, and the off diagonal blocks having entries 0 or 1, counting formulas 
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of the type q can be obtained by studying ratios of characteristic polynomials averaged over 
the classical groups. This is done in Section |21 

The topic of Section 0 is the relationship between enumeration formulas o and cca and 
certain classes of non-intersecting lattice paths, or equivalently certain configurations of vicious 
random walkers. This is motivated by the graphical representation of the RSK correspondence 
in terms of non-intersecting lattice paths nans- In Section 0 we describe connections between 
moments of characteristic polynomials and enumeration of certain classes of plane partitions. 
Finally, in Section El we revisit the question of the combinatorial meaning of the expected value 
of characteristic polynomials of random complex Hermitian Wigner matrices, and take up a 
similar study in relation to Wishart matrices. 


2 Further symmetrizations of the square 


The RSK correspondence gives a bijection between weighted nxn non-negative integer matrices, 
entries Xij weighted by (oti/3j) Xi i , and pairs of weighted semi-standard tableaux of content n. In 
the latter one member of the pair is weighted by 


a 


#l's #2's 


1 


Ctn 


• • • a 


#n's 


( 2 . 1 ) 


while the other member is weighted by GU but with the a’s replaced by the /3’s (the exponents 
in (12.11) are determined by the entries of the matrix). Summing (12.1 1) over all allowed semi¬ 
standard tableaux of a given shape n gives the combinatorial definition of the Schur polynomial, 
and this way one obtains for the generating function of the weighted matrices the identity 

T=fri 7\ XT = y ] S K {a\^ . . . , Q! rl )s K (/ll, ..., f3 n ). (2-2) 

riij=i(i - oiPj) “ 

This is well known in the theory of the Schur polynomial, and is called the Cauchy formula (see 
e.g. m ). We remark that the case of m x n matrices (m < n say for definiteness) follows by 
simply setting a m +i = • • • = a n = 0. 

The fact that the entries x^j are weighted by ( ai/3j) Xij tells us the coefficient of 

< ■■■ a **# 1 ■■■/?£* 


in (E2) counts the number of matrices for which the sum of the elements in row i equals rp, 
while the sum of the elements in row j equals pj. Because (12.21) is symmetric in the ctj’s and 
the /3/s, without loss of generality we can restrict attention to the case that r/i > ■ ■ ■ > ij n , 
Pi > • • • > Pn and thus p and p form partitions, to be denoted p and p respectively say. The 
task is then to extract the coefficient of a^/3^ (in an obvious multivariable shorthand notation). 

Only terms on the RHS of Q with |k| = \p\ = \p\ can contribute to the coefficient of 
Because of this, we can use a result of Baik and Rains |2j expressing the RHS of (EH 
with the largest part of ft restricted to be no greater than N as an average over U(N), 

n 

a*(ai,...,a n )s«(/3i, ■■■,&») = ®ueU(N) II det(/jv + otjU) det(/;v + PjU). (2-3) 

k:k,i<N j =1 
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In particular with /i, p, specified as in ed, ei agrees with El at the required order provided 


l l 

N > M = ^2jaj = ^2jbj. (2.4) 

j = 1 i=i 

The extraction of the sought coefficient from El is immediate, and we reclaim o, provided 
IV obeys the inequality El- 

We remark in passing, that since we have nm 


E t/ec/(A0 n det(/ivT OLjU') 11 det(/jvT C/ ftj') 

i= 1 i =1 


1 

(A... AO* 


(®1) • ■ • j Q!m) Ol > • • • j Ai)j (2-5) 


(in fact this equation will be of independent interest below) where IV n denotes the partition 
with n parts all equal to IV, it is also true that is equal to the coefficient of in the 

Schur function <SW«({a!j,provided N is large enough as required by (12.41) . 

From the present perspective, to extend the enumeration formula to classes of sym¬ 
metrized non-negative matrices we require a formula analogous to El for the corresponding 
partial generating function. Indeed such a formula (due to Littlewood |21| ) is known for sym¬ 
metric matrices with all elements zero on the diagonal, 

E s A (aq, • • • a) = IT ----, (2.6) 

A'even l<T<j<n 1 “ 


and for symmetric matrices with all elements even on the diagonal, 

E a *(ai, ••-«-)= n ( 2 - 7 ) 

A even l<i<j<n * ^ 

Combining El with 


n 

^MeO(N) JJdet^Ar + atjM) = ^ s A (ai,..., a n ) (2.8) 

3 =1 1(\)<N 

X' even 

and El with 

n 

® MeUSp(2N ) n det(/ 2 jv + otjM) — 'y ' s A (ai,..., a n ) (2-9) 

j=l \\<2N 

A even 

allows us to reclaim El and El respectively. The finitizations El and El of Littlewood’s 
formula are due to Rains Ill- 

Formulas analogous to El are also known in the case of symmetric matrices with no special 
constraint on the diagonal, matrices symmetric about both the diagonal and anti-diagonal, and 
matrices with a point reflection symmetry about the centre. Let us first consider symmetric 
matrices with no special constraint on the diagonal. 

For a symmetric matrix, the row sum equals the column sum so we can take for the weights 
a iPj = yjQiQj and extract the coefficient of q 11 to read off the number of non-negative integer 
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matrices with row sums equal to /i. Since each off diagonal element in position i j of a symmetric 
matrix is identical to the element in position ji, we can equivalently restrict attention to the 
lower triangular portion i < j of the matrix, and weight the elements in strictly lower triangular 
positions i < j by q%q^ and the diagonal elements by q{. The generating function for weighted 
symmetric matrices is thus seen to be 


1 




( 2 . 10 ) 


and we seek the coefficient of q^ in this expression. According to an identity of Littlewood EU 
we have that is equal to 

^2s K (qi,...,q n ). ( 2 . 11 ) 

K, 

It has been shown by Baik and Rains (2] that constrained by the size of the largest part 

of K\ can be expressed as a random matrix average, 


, • • •, Qn ) — ^ugO(n) det (/tv + U) JJ det (/tv + qjU). 

k:ki<N j =1 


( 2 . 12 ) 


With knowledge of the above results, we can now easily obtain an enumeration formula for 
symmetric non-negative integer matrices with prescribed row sums. 

Proposition 1. Consider the set of all symmetric n x n non-negative integer matrices. Label 
the row sums (or equivalently column sums) by the partition 


fl — (/ii, . . • , Hn ) — (1 1 


r). 


For N > |/i| the total number N(f of such matrices is given as a random matrix average by 


= E M60 (iV) det (/tv + M) ^(Sc^M)) 6 

3 = 1 


(2.13) 


(cf. E~3)). 


Proof. Only terms in (12.111) with |k| = \p.\ can contribute to the coefficient of q^ in (12.101) . so 
for N > \/i\ the latter is the same as the coefficient of q^ in (12.121) . Extracting the coefficient 
gives (EH- □ 

We remark that since det (/tv + M) = 0 for M G 0~(N) (because A = —1 is an eigenvalue) 
we can restrict the matrices M in (12.131) to M G 0 + (N) (the matrices U in (12.121) can similarly 
be restricted). 

Next, we note that the results of Baik and Rains in fact allow the diagonal sum to be 
prescribed. The relevant generating function is f2| 


1 


n?=i( i -«?») Ui<j0- - <Mj) 


U^y-^ SK (q u ...,q n ) 


(2.14) 
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(2.15) 


and for n\ < N this sum can be written as the random matrix average 

n 

E r/GO(Ar) det(I N + all) det(I N + qjU). 

3 = 1 

Hence, with the diagonal sum prescribed to be equal to p say and the total number of symmetric 
non-negative matrices now denoted N® p , the enumeration formula (EH should be modified to 
read 

i 

N% p = E Me0{N) Sc P (M)) JJ(S Cj -(M))^. (2.16) 

i=i 

We now turn our attention to the task of enumerating according to row sums 2 n x 2 n non¬ 
negative integer matrices [®i,j]ij=i,...,2n with a reflection symmetry about the diagonal, Xij = Xjj 
( i > j), and about the anti-diagonal, Xij = Xi^n+i-j if > 2n + 1 — j). The generating function 
for weighted symmetric matrices of this type is 


n 


(1 +Xi Qi) 


1 


(2.17) 


f = i (! - Xoft) n^(! - Wj) 

where xo = 1 (xi = 1) hi the case that the elements on the diagonal (anti-diagonal) are 
unrestricted, while \o = 0 (xi = 0) in the case the elements on the diagonal are restricted to 
be zero (anti-diagonal are restricted to be even). We know from [2] that the random matrix 
average 


Ec/ec/(jv) ^fe^H) n det(/Ar+r * t/ ) det ( / ' ?v + ( h u ) ( 2 - is ) 

correctly reproduces (I2T7D up to and including terms 0{q 2 - N+1 ) in the q y (In the case Xi = 1) 
(I2TH1) is to be interpreted as the limiting value for x\ 1 •) The following enumeration result 
is now evident. 


Proposition 2. Consider the set of all symmetric 2n x 2 n non-negative integer matrices which 
have the further constraint of being symmetric about the anti-diagonal. Introduce possible con¬ 
straints on the diagonal and anti-diagonal elements according to the values of xo and xi noted 
below With the row sums labelled by fi, we have that for N > ^\p\ the number of matrices 

in the set, N® say, is given by the coefficient of q M in the random matrix average KM- 

As our next example of a symmetry constraint, we turn our attention to the case of 2 n x 2 n 
matrices invariant with respect to reflections in the point (n +1/2, n + 1/2) (here we are thinking 
of the matrix as labelled by a grid of 2 n x 2 n lattice points : 1 < i , j < 2n}). For such 

matrices Xij = X 2 n +i-i 2 n+i-j ■ The generating function for this class of matrices is 


1 

n;w(i^) 2 ' 


(2.19) 


It follows from (12.21) and (12.51) that up and including terms 0{qf) (12.191) is equal to 


( E t/ei7(A9 n det(I jv + QjU) det (In + 

j =i 


2 


( 2 . 20 ) 
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We thus have the following enumeration result. 

Proposition 3. Consider the set of all 2 n x 2 n non-negative integer matrices with the point 
reflection symmetry Xij = X2 n +i-i,2n+i-j ■ With the row sums labelled by p, we have that for 
N > \p\, the number of matrices in the set N p say is given by the coefficient of q p in the random 
matrix average I 2.2(A) . 

As another extension of q we consider not a symmetry constraint on the matrix, but 
rather a restriction on the entries of the matrix. These we take to be either 0 or 1. For m x n 
matrices of this type, weighted by (ai/3j) Xi fl the generating function is 

m n 

n iia+<*&>■ (2 - 2i) 

i= 1 j =1 

According to the dual Cauchy identity, 

n 

JJ (l + aiPj) = ^s M /(ai,...,a ri )s At (/?i,...,/3 n ), (2.22) 

i,j= 1 V 

where p! denotes the partition conjugate to p (see e.g. EH). Analogous to (EH, we have that 
this sum restricted by the size of the largest part of p can be written as a random matrix average 

v( a i’ • • ■i a n)s il {(3 1 , ...,Ai) = E U£U{N ) JJ , (,J N 3 J-qI • (2.23) 

WUl<N j =1 QC ^ N W > 

Because of the occurrence of the dual partition in (I T22h . pi cannot exceed n for a non-zero 
contribution and so (EZI1) is independent of N for N > n. 

We want to extract from this the coefficient of a^/3^. For this purpose we introduce inverse 
secular coefficients Rc P (U) according to 

1 OO 

-—77 - —t = 'S'' x p Rc p (U). 

det (I N -xU) ^ ^ 

In terms of the eigenvalues e l>: ‘ (j = 1, ..., N) of U we have that 

Rc p (U) = h p (e ie \...,e ieN ) 

where h p denotes the pth complete symmetric function (see mi) Making use too of m we 
obtain for the sought enumeration the following result. 

Proposition 4. Consider the set of all mxn matrices in which the entries take on the values 0 or 
1. Prescribe the row sums and column sums by the partitions p and p respectively, which are to be 
written in terms of the frequency of their parts by E2P. For N > min (m, n, Yl j= i J a ji Y!j= i J b i) 
the total number of such matrices is given in terms of a random matrix average by 

i 

N°J = E ueu { N) n(Sc flU)r(R Cj (U)p. (2.24) 

3 = 1 

We remark that for to be non-zero we require p\ < m, p\ < n and \p\ = \p\. 
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3 Block matrix structures and ratios of characteristic 


polynomials 


We have seen in the previous section how generating functions for non-negative integer matrices, 
and 0-1 matrices, are related to Schur function identities which in turn are related to averages 
over the classical groups. Results obtained in 0 m tell us this strategy can also be carried 
through for classes of block matrices which relate to averages over the unitary, symplectic and 
orthogonal groups. 

In relation to the unitary group, let the block structured {k\ + l\) x [k 2 + h) non-negative 
integer matrix 


A B 
C D 


(3.1) 


be such that A is a k± x k 2 non-negative integer matrix; D is a l\ x l 2 non-negative integer 


X = [ Xij \ are weighted by gA J with 


9ij — < 


i a l\ 

xk 2 0-1 matrix. 

Suppose the entries of the block 

a ilfj 

for X = A 


Pi$j 

for X = D 

(3.2) 

oiidj 

for X = B 

Pm 

for X = C. 



Such matrices have for their generating function 


fcl k2 

nn 


t=1 3= 1 


1 h h 

1 - n 

j i=i 7=1 


1 

1 - PiSj 


nnd+^onfid+AT,), 


i=lj=l 


i= 1 j=l 


thus containing as special cases both the LHS of D and (12.211) . The coefficient of aT/3^ 'fb v 
tells us the number of matrices HU with prescribed row and column sums. Moreover general¬ 
izations of both the Cauchy identities (IE3> , dZZ> , and their finitizations n . are known 

00, allowing for the counting function to be expressed as a random matrix average. 

To state these generalizations requires introducing the functions m 


(c^l Oi-k )/^1 j • • • fil) d.6t — l)l<i,jf<Z(A) 5 (3.3) 

where a& denotes the coefficient of x k in 

U l j= i(! + ^-a) 

nLii 1 -) 

The HSa are referred to as the hook Schur functions. In terms of this the generalization of the 
Cauchy identities is 


^ ) RSa(qu,..., cxki: Pi 1 • ■ ■ Pi\ )HSa( yi,..., 7 a ?2 , ^1, • • • ^2) 
A 


k 1 &2 


fcl I 2 


h k2 


nn T v^nn T v M nn( i +“-«nn( i +^) 

J i =1 7 = 1 J i=l 7 = 1 


i= 11=1 


*=11=1 


(3.4) 










while the generalizations of their finitizations is 


E rii=i det (In + atjU) n -li det (I N + 7 jU) 

U£U{N} nlu det(/jv - PmU) nti det (/tv - 

^ ^ (cTl , . . . , CK/ci 5 /^1 1 * ' ' )HS_\ (^1 j • • ■ ; r //C2 ; ? ' ' ■ ^2 )* (^*^) 

M<N 

Arguing as in the derivation of the counting formula o given in the first three paragraphs 
of Section ID the following generalization of (II .‘ 2 jl and (1A241) is immediate. 

Proposition 5. Let a = (ai,...,a fel ) , b = ( 61 ,..., b h ), c = (ci,...,c fc2 ) , d = (d 1 ,... ,d h ), 
where the a,j, bj, Cj, dj are nonnegative integers. From these arrays form partitions 

a* = U ai ---C 1 ), A = (1 61 ^ = (1 C1 ---^ 2 ), i> = (i d '---l d 2 h ). 

Let denote the number of matrices 13. 1\) with 

row (A, B) = fj,, row(C, D) = v, col (A,C) = ju, col (B,D) = D 


where the notation row(X, Y ) refers to the row sums across X and Y given they are horizontal 
neighbors in a block matrix, and col(X, Y) refers to the column sums down X and Y given they 
vertical neighbors in a block matrix. 

For 

( k\ l\ k,2 I 2 

jbj^jacj^jdj, 

111 1 


we have 


k\ li /C2 f-2 

®u N n(Sc,(M))- n (Scj(M)) bj (Rc m (M)) Cm J] (Rc n (M)) dn = N^. (3.6) 

i= 1 j =1 m= 1 7i—l 


As reviewed in the Introduction, a relationship with averages over the orthogonal and sym- 
plectic groups comes about when the matrix is constrained to be symmetric. Thus in ED we 
must take 

ki = k 2 = k, l 1 = l 2 = l, A = A t , D = D t , C = B t . (3.7) 

Correspondingly, for i j, the weights EH) are to be replaced by 


yfajaf for X = A 

for X = D 

for X = B 

\j PiOLj for X = C. 


(3.8) 


As in the meaning of CD and CD, the constraint on the diagonal elements will determine 
whether the relationship is with an average over the symplectic group, or an average over the 
orthogonal group. It turns out that the constraint relevant to the symplectic (orthogonal) group 
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is that all elements on the diagonal of A be even (zero), while those on D be zero (even). The 
corresponding generating functions are therefore given by 


k l 


gSp «“-}-^» = n n T^nn™.,) 

J i =1 j=i 


\<i<j<k 3 1 <i<j<l 


3= 
k l 


G°(w.{ ft }) = n n iz^nno+o*) 

^ 1<2<7</ ^ 2=1 .7=2 


1<*<?<2 


2=1 J =2 


In relation to these generating functions, one has as generalizations of the Littlewood iden¬ 
tities EH) and (12771) gj 


E HSa(q:i, ...atift,... ,Pi) = G Sp ({ai},{/3j}) 

\ even 

E HSa(c>!i, ... a k -,/3i, ■ ■ ■ ,A) = G°({«i}, (d})- 


A' even 


And with the largest part of the partitions in the sum restricted, the LHS’s have the finitizations 


E. 


nil det(/ 2 AT + a-iU) 

u £U Sp( 2 N) pw—TT77-T77T 

n ;= idet(/ 2 iv - PjU) 


E, 


UzO(N) 


nti det(Ijy + ajU) 
n'=, det(l v - PjU) 


E HSa(q;i, • • • CKfc;/3i, ■ ■ ■ ,A) ( 3 -9) 

Ai<2jV 
A even 

E HSa(o;i, • • • a k ;/3i,... ,A). (3.10) 

l(X)<N 
X' even 


As a consequence, the following enumeration results hold. 

Proposition 6. Let a = (oi,..., a*,) , b = (bi, ..., bi), where dj, bj are non-negative integers, 
and from these arrays form partitions 


H = (l ai • • • k ak ), u = {l bl ■ ■ ■ l bl ) 


Consider block matrices El with constraints and the further constraint that all diagonal 

entries of A are even (zero) while all those of D are zero (even). Let N^f) (AT®,) denote the 
number of such matrices with 


row(A, B) = col(A, C) = n, row (C, D) = col(i?, D) = v. 


We have 


k i 

pv ®'M£[/Sp(2A r ) n< sc i< M ))“‘ n<**<*»‘ 

i=l j= 1 


N pu — E MeO(7V) 


k l 

n(sc,(M))-n(Rc,(M)).-. 

3 =1 


i=1 
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Figure 1: Example of non-intersecting lattice paths corresponding to returning lock step 
vicious walkers. For the first (last) n steps the walkers must move to the right (left) or 
stay stationary. The weight of the configuration shown here is a\a\a\a^(3l faPzPt- 

4 Relationship to non-intersecting paths 

The purpose of this section is to relate the random unitary matrix average in El and the 
unitary symplectic average in El to configurations of weighted non-intersecting lattice paths. 
This then allows us to give combinatorial interpretations to El and El relating to non¬ 
intersecting lattice paths rather than integer matrices. That such interpretation are possible 
can be anticipated from the RSK correspondence: the semi-standard tableaux therein have a 
well known interpretation in terms of weighted non-intersecting lattice paths (see e.g. ffi)- 
However neither in the case of El nor El will our non-intersecting paths correspond to the 
conventional ones related to semi-standard tableaux, and so a separate discussion is warranted. 

In relation to El, mark points on the x-axis at x = 1,... ,N. Move each point to the line 
y = 1 according to the rule that each x coordinate must either stay the same (weight unity) 
or increase by one (weight au), with the proviso that all x coordinates must remain distinct. 
Connect the points between y = 0 and y = 1 by segments, which must either be vertical (weight 
unity), or right diagonal (weight ai). Repeat this procedure a total of n times, with each right 
diagonal segment at step j weighted by atj. 

After step n, perform another n steps, but now with the segments either vertical (weight 
unity) or left diagonal (weight fon+i-j hi step n + j). The segments again must not intersect, 
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and are further conditioned to return along y = 2n to the same x coordinates (x = 1,... ,N) 
as they began (see Figure Q for an example). The resulting non-intersecting lattice paths are 
equivalent to a special case of the lock-step model of vicious random walkers mm- For general 
initial positions along y = 0 (4°\ ■ ■ ■ , 1$ say), and final positions along y = 2n (l±,... ,In say), 
the generating function G 2n for the weighted paths can be written as an N x N determinant 
according to 


G 2n (l { °\...,l ( £ ) -J 1 ,...,l N ) = det 


92n{l? ) -,h) 


J j,k=l,...,N 


(4.1) 


where 

<?2n(Z (O) ;0 = ^~ r n(l + a J -e- < ®0(l + /3 7 -e w Oe“ i(, “' (0)) ^- (4.2) 

71 J - n j =l 

From the well known identity 


det 


- 1 
.2tt 



h{6)e-^- k)e dd 


j,k=l,...,N 


(nw)^, 

3 = 1 


(4.3) 


we see that in the case of interest (initial positions = final positions, all one unit apart), we have 


G 2 n(q 


( 0 ) 


/(°).K 

> > *' 1 > 


Jn) 


] 3 

0 = 1 , 


— ^ugu(n) n det (In + ctjU) det (In + PjU). 

3 = 1 


The following interpretation of the average over U(N) in (11.21) is now evident. 


(4.4) 


Proposition 7. Consider the set of all non-intersecting paths of the type depicted in Figure 
Q Impose the further constraint that the number of right diagonal segments between y = j — 1 
and y = j (j = 1,... ,n) is equal to fij, and the total number of left diagonal segments between 
y = 2n — j and y = 2n — j + 1 (j = 1,... ,n) is equal to fij. The total number of such lattice 
paths is given by the RHS of m. 

We note that in the setting of Proposition 0 at most |/r| = (which must equal \fi\) 

different walkers, counted to the left from x = N, can move. Thus we see immediately that o 
must be independent of N for N >\p\. 

Let us now turn our attention to the non-intersecting lattice path interpretation of (EH). For 
this mark points on the x-axis at x = 1,..., N. Also count steps j = 1,..., 2n as the points are 
moved to the lines y = 1 ,..., 2 n in order, and the corresponding segments are drawn to connect 
the points. Let this process proceed by the rule that the lattice paths must not intersect, and 
that at odd numbered steps the segments must be either vertical or right diagonal, while at 
even numbered steps the segments must either be vertical or left diagonal. Weight the diagonal 
segments in steps 2 j — 1,2 j by dj,bj respectively. It is further required that all points remain 
to the right of the y-axis, which can be thought of as a wall (see Figure El f° r an example). 
For general initial positions along y = 0 (l^\...,lff) and finishing positions along y = 2 n 
(l \,..., In), all to the right of the wall, the generating function G ^ 11 for the weighted paths has 
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Figure 2: Example of non-intersecting lattice paths with initial (final) spacings along 
y = 0 (y = 2n) one unit apart starting at x — 1. All paths are restricted to x > 1, and 
segments can be vertical or right diagonal for odd steps, and vertical or left diagonal for 
even steps. The weight of the configuration shown here is afafalal. 

a determinant form m ■ In the special case that {a*} = {bi} = {a*} (in any order), and that the 
initial and final positions are spaced one unit apart starting at x = 1, this determinant can be 
expressed as the random matrix average m nn, thus providing the following interpretation 
of the latter in terms of non-intersecting paths. 

Proposition 8 . Consider the set of all non-intersecting lattice paths of the type depicted in 
Figured Impose the further constraint that the sum of the number of right diagonal segments at 
step 2 j — 1 and the number of left diagonal segments at step 2 j is equal to yj. The total number 
of such paths is given by the RHS of db 

5 Relationship to plane partitions 

In this section the random unitary matrix average (12.51) . the random unitary symplectic average 
(HH), and the random orthogonal average (12.151) will be related to the enumeration of certain 
classes of plane partitions muni- a plane partition V is a finite set of lattice points {(i, j, k )} C 
N 3 with the property that if (a, b,c) £? and 1 < i < a, 1 < j < b, 1 < k < c, then ( i,j , k) £ V. 
A plane partition is symmetric if ( i,j,k ) £ V if and only if ( j,i,k ) £ V. The height of stack 
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(i,j) is the largest value of k for which there exists a point ( i,j, k) in the plane partition. 

The study of plane partitions was initiated by MacMahon m who proved that the generating 
function for plane partitions fitting in the box 


B(a, b , c) = lfe)|l < i < a, 1 < j < b, 1 < k < c} 


(5.1) 


is given by 


a b c 


nnn 


1 _ qi+j+k -1 

X _ qi+j+k- 2' 


(5.2) 


Taking the limit q —> 1 gives for the total number of plane partitions fitting inside B(a,b,c), 
#V(a,b,c) say, the evaluation 


#V{a,b, c) 


a b c 


=nnn 


i + j + k — 1 
* + j + k — 2 ’ 


(5.3) 


We can express the generating function (fT2l) for plane partitions in terms of Schur functions. 
Thus from the combinatorial definition of Schur functions sa(xi,..., x n ) as a sum over weighted 
semi-standard tableaux of shape A and content {1,... , n}, it follows that 


(a+c a+c—1 
s b a [Q , Q , 


■ ■,q) 


where b a denotes a partition with a parts all of which are equal to b, is the generating function 
for plane partitions strictly decreasing down columns with exactly a rows each of length b and 
with the largest stack of height less than or equal to a + c. The strictly decreasing constraint 
can be eased by removing a — i + 1 from the boxes in row i to bijectively obtain a plane partition 
which is a subset of B(a,b,c). Consequently an alternative expression for (15.21) is given by 

q -ba(a+l)/2 Sba ( q a+c^ g a+c- 1, _ _ ^ (5.4) 

Making use of (12.51) shows that for a = n, c = m and b = N, (IP1) can be expressed as 

m n 

^ueu(N) n det(ijv T OijU ) 11 det(/jv T U/ftj') 

j =i j =i 

with ati = q n+l and (3i = q l . Taking the limit q —> 1 gives the sought relationship between (15.31) 
and a random matrix average. 

Proposition 9. Denote by P^{z) the (reciprocal) characteristic polynomial of M E U(N), so 
that Pm{z ) = det (Ijy — zM). For \z\ = 1 we have 

#V(a,b,c) =E t/G[/(c) (p A * / ( 2 )) 6 (jVi)) a . (5.5) 

We remark that in mi the case a = b of the RHS of Q has been given a combinatorial 
interpretation involving two-rowed lexicographic arrays. 
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The random matrix average is an example of a class of multi-dimensional integrals 

possessing gamma function evaluations (see e.g. m)- This implies 


#V(a,b,c) 


tt r(a + b + 1 + j)T(2 + j ) 
i} Q T(a + l+j)T(b + l+j) 

G(l + a + b + c) G(l + a) G(l + b) . . 

G(1 + a + b) G{l + a + c)G(l + b + c) [ ' 


(5.6) 


where G{z) is the Barnes G-function, related to the gamma function by the functional equation 
G(z+ 1) = T(z)G( 2 :), G(l) = 1. In contrast to the formula (15.511 . the latter formula is well suited 
for asymptotic analysis. 

Consider now symmetric plane partitions fitting inside the box do with b = a. The 
generating function for such plane partitions, with the additional constraint that the heights of 
all stacks on the diagonal are even and bounded by 2 c, has the product form IS 


n 

l<i<7<a 


1 _ q i +J + 2c 

1 - q i+ J 


(5.7) 


Here only points on and above the diagonal are weighted. Denoting their total number by 
#Ve ym (a, 2 c), we see by taking the limit q —■> 1 that [7 

#Pf m (a,2c) = n Z + 1 +2C - (5-8) 

i Aj< a 1+3 

The generating function Q and thus counting formula (15.81) can be expressed in terms of 
Schur functions by making use of a bijection between tableaux and symmetrical plane partitions. 
Consider then semi-standard tableaux of content {l,...,a}. Suppose furthermore that each 
row is of even length, and the first row is constrained to be less than or equal to 2c. From 
this construct the diagonal and upper triangular portion of a symmetrical plane partition by 
associating with grid points i < j, stacks of height 


hi,j = #^5 + #(* + l/s + • • • + #(a + i — j)'s in row i of the tableaux. 

By weighting each square labelled j in the tableaux by q J , we see that each stack at grid point 
(i,j) is weighted q hi ’G and furthermore on the diagonal this weight is where Aj is the length of 
row i of the tableau. It follows that (E3> can be expressed in terms of Schur functions according 
to 

' s A(<f,g ,a- \ ■■•,<?)• (5.9) 

AC(2c) a 
A even 

Recalling gives the sought relationship with an average over the unitary symplectic group. 

Proposition 10. Denote by P^f(z) the characteristic polynomial of M E USp(2c ) so that 
pSf(z) = det(A/ 2 C — M). We have 

#Pf m (a,2c) = E Met/5p(2c) (pg > (-1))“. 
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Analogous to ra . the above average is of a type which admits a gamma function evaluation 
(see e.g. na) This allows (EEJ) to be written 


#Pf m (a,2c) 


2 2 ca TT r (l + C +j)r(^ + a + j) 

j=i r(i + c + a + j)r(i + j) 

o2 C a ^( 2 + 2c ) G(2 + c + a) (2(§ + a + c) <?(§) 
G(2 + c) G(2 + 2c + a) G(f + a) G(|+c) 


Finally, consider the same weighted symmetric plane partitions as described in the paragraph 
below ra . but with stacks now restricted to heights < c and without the constraint that the 
stacks on the diagonal be even. The generating function is 


s A(<f,g a 

AC c a 


In an obvious notation, after taking the limit g —>■ 1 and recalling dm we obtain 

W sy m (a,c) = E MGO(c) (p°/(-l)) a+1 . 


6 Counting formulas associated with the characteris¬ 
tic polynomials of random Hermitian and Wishart 
matrices 

Let X be an n x n Hermitian matrix. Let the diagonal elements xa be chosen independently 
according to a probability distribution V\, with the property that 

E Vl xu = 0 . ( 6 . 1 ) 

Let the upper triangular elements Xij , i < j, (which may be complex) be chosen independently 
according to a probability distribution V 2, with the properties that 

^v 2 Xij = 0, Ex, 2 \xij\ 2 = of. (6.2) 

The following result, generalizing a number of results presented in [5] 1 , gives the expected value 
of the characteristic polynomial of X in terms of the classical Hermite polynomial H n (x). 

Proposition 11. Let X be Hermitian and specified in terms ofT>\ and T> 2 as above. With the 
monic rescaled Hermite polynomial specified in terms of the classical Hermite polynomial by 

h n (x) = 2-"' 2 H n (fL') 

we have 

Evi_,v 2 det(A/jv — X) = cx^/t/v (—)• (6-3) 

A CT'2 ‘ 

1 We take this opportunity to correct the statement of Theorem 15 in [3]. It should read: 

E hn(Pm{ x )) — h^\x), where hffl are orthogonal polynomials with respect to the weight (t — x) 2k e~ t2 . 
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Proof. By definition 


N 


det(XlN — X) — ^ £(P)Y\.(^i,p(i) ~ Xi,p(i)), 
PeS N 1=1 


(6.4) 


where e{P) denotes the parity of P and 



The specifications dSJ) , (tO) tell us that the only non-zero terms in (IQl) after averaging 
over T>i,T >2 are those for which P consists entirely of fixed points ( P(j ) = j) and 2-cycles 
(P(j 1 ) = ,]2 and P(j 2 ) = ji, j 1 7 ^ jfi)- Let there then be N — 2 j fixed points and j 2-cycles. 
Such permutations have parity (— l) 3 . Each fixed point contributes a factor A, while each 2-cycle 
contributes a 2 . As the number of ways of choosing the (N — 2 j) fixed points and the j 2-cycles 

'N\ (2 j)\ 

■ -j) 2Jj! 

we see that 


E© 


i^ 2 det(A I N -X) = E(-l ) j (^)^ AAr - 2 V 2 2 T 


But this is precisely the power series expansion of the polynomial in question. 


□ 


Next we turn our attention to the mean characteristic polynomial of so called chiral matrices, 
that is matrices of the form 


K : = 


0 rax 7i 


X , 
0 ri 


nxp 


(6.5) 

Jpxn Vpxp 

where we require that n > p. These matrices have exactly n — p zero eigenvalues, with the 
remaining 2 p eigenvalues given by ± the positive square roots of the eigenvalues of the non¬ 
negative matrix X^X. We specify that the elements of X be chosen independently according to 
a probability distribution V with the properties 


E T>Xij = 0, Ex>\xij\ 2 = a 2 . 

The analogue of Propositional can readily be deduced. 


( 6 . 6 ) 


Proposition 12. Let K be a random chiral matrix as specified above. We have 

E v det(A I n+P - K) = pla 2 P\ n -PL;-P((\/a) 2 ), (6.7) 

where L^(x) denotes the classical Laguerre polynomial. 

Proof. Because of the first specification in (IQll . we see that in the analogue of (16.41) for 
det(A/ n _(_p — K ), the only terms after averaging will again result entirely from fixed points and 
2-cycles. The fact that K has n — p zero eigenvalues implies we require there be a minimum of 
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n — p fixed points. Thus we must consider the cases that the number of fixed points is equal to 
n + p — 2 j and the number of 2-cycles is equal to j for each j = 0 ,p. 

Not all permutations with these specifications give a non-zero contribution. For the latter, 
because of the zero blocks in (EH), in relation to the 2-cycles (.7 1 ,7 2 ) we require 

3i e 77.}, j 2 G {n + 1,..., n + p}. 


Hence we must choose ( n — p) + (p — j) fixed points from the first of these sets and (p — j) from 
the second. The number of distinct ways to do this is 


n 


P 


-n-j; \p — j > 

The j 2-cycles can then be chosen in jl different ways. Since again each fixed point contributes 
a factor A, while each 2-cycle contributes a factor a 2 , we see that 

p 


Epdet(A In+p-K) = a 2 p X n ~ P Y^ 

j =0 

p 

= P \a 2p X n ~ p J_l 
j =0 


n 


P 


n- 3 J \p- 3 


X\ 2 (p~j) 


Mi) 


n 


p- jJ j\ \ (J 


1 f x y j 


The sum in this expression is precisely Lp p ((X/a) 2 ). □ 

Due to the relationship between the chiral matrices X^X as noted below EH- we have the 
following result for the expected value of the characteristic polynomial for Wishart matrices. 


Corollary 1. Let X be a n x p matrix with elements independently chosen according to the 
distribution D with the properties \ 6 .b\) . We have 

E-p det(A I p - X ] X) = pla p L™- p (X/a). (6.8) 
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